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We develop a renormalization group (RG) procedure that includes important system-specific 
features. The key ingredient is to systematize the coarse graining procedure that generates the RG 
flow. The coarse graining technology comes from control and operator theoretic model reduction. 
The resulting "generalized" RG is a consistent generalization of the Wilsonian RG. We derive the 
form of the projection operator from the dynamics of a nonlinear wave equation and renormalize 
the distribution of initial conditions. The probability density of the initial conditions is chosen to 
be the Boltzmann weight for a standard 0^-theory. In our calculation, we find that in contrast to 
conventional implementations of the RG, naive power counting breaks down. The RG-equations 
that we derive are different from those derived from the conventional RG. 

PACS numbers: 

I. INTRODUCTION 

A recurring theme in the study of complex and biological systems is that the systems of interest are open. The 
nature of the disturbances rendering the system open may be stochastic, structured, or both. Additionally, the 
stochastic nature of chemical networks, especially those with molecular-types in low abundance, ensures that any 
intrinsic description must account for the structure and type of uncertainty. Intuition from closed systems, with 
the possible exception of glassy systems |20], suggests that detailed questions are often too difficult to answer and 
that asymptotic questions are more amenable to theoretical treatments. While for many open systems knowledge of 
the temporal asymptotics is both interesting and important, for others such knowledge may not carry any content 
or even exist. However, just because temporal asymptotics may not exist for open systems, this does not imply 
that there aren't robustly discernible features of open systems. It simply means that the appropriate coarse-grained 
description(s) must be informed by the structure of the external and internal disturbances and uncertainties. 

A large class of open systems may be represented by equations of the form 

^^ = /(M)-H«(x,o, (1) 

where </? and u are vectors/functions in possibly infinite dimensional spaces. Without u, equation ([T]), represents 
a closed system or a model that describes a system with perfect certainty. With u, the system is inherently open 
subject to additive uncertainty, u represents generic driving as well as a possible noise source to which the system is 
exposed. If u is generated by a continuous stochastic process, this imposes a specific structure on the noise. Similarly, 
constraining u to belong to a particular function space but remaining otherwise arbitrary also imposes a specific 
structure on the noise. The system, as described, may be mapped to a field theory via generalized Martin-Siggia- 
Rose (MSR) /closed-time-path (CTP) methods [^H, [ifl, H^- In particular, the path-integral representation for the 
probability for the system to be in state at time t given that the initial state was 0o at time to is given by [1^ 

P((/.,#o,to) = j ^V'expl-^^ dTd^f]5[u^dr'f + f{ip) + (po5{T^to)], (2) 

When the initial state is only known probabilistically, then the probability for the system to be in state at time t 
is given by 

j DcpoP{^oMP{M^oM (3) 

We introduce this representation because later we will explicitly utilize it in the context of applying the renormalization 
group (RG) to the initial conditions. It is important to note that upon combining Equations ^ and ([31), integrating 
out does nothing to dress the distribution P{(j)o,to). This simply reflects the fact that the initial conditions are 
imposed as a constraint. On the other hand, integrating out 0o would then add stochastic contributions to Equation 
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Suppose u is an arbitrary input to the system that is square-integrable in time. We consider the states or regions 
in phase space that are most accessible via driving to be responsible for describing the essential characteristics of the 
system. This is analogous to the energy landscape picture in statistical mechanics where fluctuations govern which 
states contribute the most to the statistics of the system. We use this inherently open systems perspective of the 
importance of states to specify how to coarse grain and, consequently, to generate RG equations. The key step in 
generalizing the RG lies in ascertaining how to coarse grain. 

Although it is already known that the RG is not a black box routine, the purpose of this work is to make it more 
algorithmic. It is easy to be misled into thinking that the RG already is algorithmic because its key ingredients are 
coarse graining and rescaling the system variables [H, [1, [1, Q . However, fully algorithmic implementations of the RG 
fail for large classes of problems because it is not possible to ignore the physics of a system and expect to obtain 
meaningful results. Capturing the essential physics requires isolating the appropriate models and the structure of 
perturbations and uncertainties. Consequently, this process is system-specific. Additionally, the scale on which the 
physics is observed must be specified. For instance, for bosonic theories the long-wavelength physics is investigated, 
while for fermionic system it is the physics near the Fermi surface. These considerations suggest that the primary 
obstacles to automation are the model identification and coarse graining processes. In this article we present a 
RG procedure that accounts for these system-specific obstacles and apply it to renormalize the initial conditions of a 
nonlinear wave equation. We specifically focus on a nonlinear wave equation as opposed to a reaction diffusion equation 
because the extra time derivative increases the number of Green's functions that can contribute to the dynamics of the 
system. As we will see, this can have a significant impact on coarse graining. We do not renormalize the dynamical 
equations because the structure of the dynamical action renders the perturbative loop corrections uninteresting. The 
general form of how projection operators, generating the RG, act on the dynamical fields is discussed in p6j . 

The obstacles mentioned above arise from the (mis) identification of observables. We use techniques taken from 
operator theoretic interpolation theory and control theory to systematically identify observables that respect the 
system dynamics. The conventional implementation of the renormalization group relies upon discerning equilibrium 
observables at predetermined scales. Coarse graining nonequilibrium systems must respect the nature of the dynamics 
and the structure of the uncertainty. Emerging fields like systems biology are producing many new problems that are 
begging for such a treatment [l7j . This work provides a step towards the development of a framework that will 
be applicable to such systems. 

While equation ([1]) addresses the effects of perturbations to the system, it does not account for the consideration 
that it may not be desirable to describe every observable. For instance, complete characterization of a system over 
all scales may carry significant computational overhead. Furthermore, due to experimental constraints, it may not 
be possible to measure all observables either. Consequently, it is beneficial to explicitly include the possibility of 
multiscale or constrained observation. One of the simplest cases is when the observables are linearly dependent on 
the (/3-observables. The description of such open systems takes the form 

V = fi'P) + Bu, 

= C^, (4) 

where ip reflects that only some subspace or, more generally, subset of phase space is directly measurable. The 
operators B and C respectively specify the structure of how noise may enter the system and which states are considered 
as observables. More generally the observables may be nonlinear functions or functionals of ip, however without 
mapping the problem to a yet higher dimensional description, that complication is beyond the scope of the analysis 
presented herein. Considering coarser observables "0 influences the relative importance of the internal states, if, 
because many of them simply either do not or only marginally contribute to the physical observables ip. Hence, the 
choice of observables can strongly influence how models describing a system ought to be coarse grained. For instance, 
in the original analysis of Feynman and Vernon or Caldeira and Leggett ^18j characterizing a particle in a heat 
bath (albeit in a quantum mechanical context), the microscopic degrees of freedom contributing to the heat bath were 
systematically projected away thereby leaving their effective influence on the particle. Implicitly, in their analysis a 
choice of C was made and consequently led to the derivation of effective stochastic equations describing the physics 
of the chosen observables. 

The outline of the paper is as follows. In Section [H] we introduce the coarse graining procedure. We use the 
procedure in Section IIIII to identify the projection operator for diffusion and wave equations. In Section IIVI the 
projection operators are used to generate RG-flows. We derive the associated RG-equations and contrast the RG- 
fiows. Our conclusions are given in Section|Vl Two Appendices provide detailed calculations related to the derivation 
of the projection operators and the RG-equations. 
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II. COARSE GRAINING VIA BALANCING 



We coarse grain by retaining the states contributing most to the dynamical response of the system. While techni cally 
different, this work is in same spirit as that by Chen et al [1^ and more recently by Degenhard and Rodriguez [26l.l27l|. 
As a starting point, we coarse grain equation ^ from its linearization about a particular solution with u = 0. The 
reasoning for this approach follows many of the same arguments why linear-response theory has been used with 
such frequency. Specifically some of those reasons are that 1) the approach is amenable to calculation, 2) it is often 
known that the system operates about a particular equilibrium, and 3) for large-amplitude noise, the linear-response 
dominates since the system then does not get trapped in local equilibria. As a simplification, we only consider 
linearizations about dynamical steady states. In systems biology, similar approximations are ubiquitously used when 
describing biochemical networks (IGj . The linearizations are generically described by 

If = A(p -\- Bu, 

^ = C^, (5) 

where B and C are the same operators as those in (j4]), = cpeq + >f where ipeq is an equilibrium solution, and 
A = f^v/lvcg- Associated with equation ^ are invariants known as Hankel singular values (HSV's) [2l|, [2^. These 
invariants have many nice properties, but their most significant being that they give explicit information about the 
Green's functions (i.e. dynamical evolution operator) for the system. These invariants are most easily calculated by 
solving for positive operators X and Y, that are determined by the equations 

^^AX + XA^ + BB^; X(0) = 0, (6) 
dtf 

dY 

A^Y + YA + C''C; r(0) = 0. (7) 



dt 



f 



If we define another operator W by 



= XY, (8) 



then the HSV are nonnegative real numbers, (Tmax > > cTmin that comprise the spectrum of the operator W. 
HSV's provide a precise measure of the error incurred by approximating the effect u has on ^p with reduced order 
models. The HSV's may be interpreted as supplying a measure of the importance of the internal states (p. If W 
is invertible, it is always possible to find a coordinate system, called balanced coordinates, such that X = Y = 
diag(tTmax, ■ • ■ J fmin)- When equation ([S]) is transformed to balanced coordinates, the best reductions are those that 
project out the states corresponding to small HSV. In other words, the ordering of the HSV, at least locally around 
an equilibrium configuration in phase space, specifies how to coarse grain a system. An in depth treatment of this 
material may be found in 22]. It is also sometimes possible to "balance" the full nonlinear system [2^ . 

The RG can easily be adapted for HSV-based coarse graining. Operator theoretic approaches to the RG [1, 0, [1| 
demonstrate that coarse graining in the Wilsonian RG is equivalent to multiplying operators or states by projection 
operators. The essence of this work is to use HSV's to identify the projection operator. As before, suppose that 
K is a vector index that orders the HSV's ctk for equation ([5]) from largest to smallest. A generalized Wilsonian 
RG procedure is obtained by: 1) transforming both the initial conditions and the dynamical variables to balanced 
coordinates, 

ifinyt) = J _R(k, x)(p(x, t)dx, (9) 
so that the distribution P{(j>, t) takes the form, 

P(0,t)= J D(boDipD<fie^p{~Sdyn (g,M,{(^},{</)o})-^(g,{0o})) 
= J D$oD^D0J exp (^^Sdyn (g, {^}, {0}, {M) - S (g, {0o}) ) , (10) 

where Sdyn is the action representing the dynamics, S is the action of the initial conditions, g is the original set of 
coupling constants/functions, (p is the dual field that arises from representing the functional-5 in equation ([2]) as an 
exponential, J7 is the Jacobian from equation Q, and g is the resulting transformed set of coupling constants; 2) 
integrating out K-shells about CTmin analogously to wave-vector shells; and 3) rescaling k and (p appropriately. An 
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interesting but technically challenging variant of this procedure is to integrate out (TK-shells instead of K-shells about 
CTmin- Integrating out a single CK-shell may entail integrating out an entire subspace in x-space because cr„ does not 
necessarily respect spatial dimension. The technical challenge lies in rescaling CTk- It is not clear that rescaling CTk will 
recover the full K-space thereby generating a meaningful RG. 

A standard interpretations of the large wave-vector cutoff in the RG is that the inverse-cutoff is proportional to 
the smallest length scale of the system 0, The idea of projecting out the large wave- vector physics down to 
the small wave- vector physics has the interpretation of homogenizing the system to its continuum, large-scale limit. 
However, interfaces that are smooth on short spatial-scales may actually appear sharp when viewed on large scales. 
Representing such effective, sharp interfaces, in terms of Fourier (wave-vector) modes then requires a huge number 
of modes (the largest mode being the cutoff). Other representations, like wavelet representations, are much better 
equipped to describe the singular nature of such interfaces. Mathematically this means that the cutoff is specifically 
chosen because it provides us with a starting point for coarse-graining and an ordering relation that determines in 
which direction to coarse grain. The cutoff can be thought of entirely in approximation/interpolation theoretic terms. 
The initiation of coarse graining starts with modes contributing the least to the phenomena of interest. The ordering- 
relation results from the ordering of modes by their contribution to the approximation. CTmin is a generalization 
of the standard wave- vector cutoff and the ordering-relation, determined by the ordering of the HSVs (a^), is the 
generalization of standard scale-ordering (small scale to large scale). The RG has been studied in the context of 
approximation/interpolation theory before. RG methods have been related to theory of Fade approximants by Baker 
PH . to wavelet theory by Battle [12, 13J, and to subspace/Lanzcos methods by White jl4 il5i|. 



III. IDENTIFICATION OF THE PROJECTION OPERATOR 



Here we apply the aforementioned coarse graining procedure to identify the projection operator for linear diffusion 
and linear and nonlinear wave equations. Since the procedure requires us to expand about a dynamical steady state, 
if we expand about the trivial solution, then the nonlinearity does not change anything. In this section we find that 
by considering slow observables for the diffusion equation and wave equation, projecting out large wave-vector shells 
is appropriate. When we additionally consider the fast, kinetic observables for the wave-equation, we arrive at quite 
a different result. This results from the extra time derivative of the wave-equation and is the very reason why we 
consider the nonlinear wave-equation in this paper. 



A. Linear diffusion equation 

We first consider the (driven) diffusion equation 

dt(l) = DV^(l) + ju. (11) 

In this example, B = C = 1, A = DV^ and we take tf — > oo. By considering a stable system over an infinite time 
horizon, we only need to solve the Lyapunov equations, 

AX + XA'' + I3B^ ^0, (12) 
A'lY + YA + C^C ^0, (13) 

instead of equations ([6])- ([7]). By taking the Fourier transform, equations (fT2|) - p^ become 

2i?|k|2Xk + |7|2 = 0, (14) 
2D\k\^Yk + l = 0. (15) 

From equation ([8]), it follows that in balanced-coordinates W is given by 

Here the index, k for the HSV's is just |k|. Thus, for the diffusion equation, the most important states are those 
that correspond to small wave-vector. Thus, local coarse graining is appropriate because the smallest observable 
"fluctuations" are due to the short-wavelength physics. The smallest error is incurred by projecting out large wave- 
vectors. We have thus identified the projection operator for the diffusion equation and it exactly coincides with what 
we would expect. Diffusive dynamics spatially homogenizes disturbances, thus, intuitively we already know that local 
coarse graining is desirable. 
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B. Linear wave equation 



As a second example, we consider the (driven) linear wave equation 



(17) 



By the units of u, it represents a true force acting on (p. This, in addition to the fact that (f) is the "measurable" 
quantity, implies that we have isolated our attention on 0-based observables. This choice is strongly influenced by 
equilibrium statistical mechanics and thermodynamics. We have completely neglected tt, the field conjugate to 0, 
that represents the kinetic contributions to the system. When posed as a set of first order equations, equation (fT7|) 
becomes 



(18) 



By smoothing out the timc-cutoff at with a damped exponential in the integral representation of the solution of 
equations ([6]) and ([7]) the problem simplifies to solving Lyapunov equations. This smoothing process is also known as 
exponential discounting. With the given form of B and C in this problem, we find that the matrix of HSV's, W , is 
approximately given by 







" 10" 






dtTT 




u^V^ 7 




vr 


y 




1 




u 



w. 



bal 



l7l 



4ati|k| 



>I. 



2x2, 



(19) 



where I2x2 is the 2x2 matrix identity, a 1/tf, and (g) is the tensor (dyadic) product. The calculation that produces 
this result is a special case of the calculation in Appendix [A) As with the diffusion equation, short-wavelength physics 
does not significantly contribute to the response, so locally coarse graining is appropriate. In the examples considered 
here, the dynamical equations that specify the projection operator are linear. Thus, no approximation has been made 
and we expect that the corresponding projection operators are globally valid in phase space. 



C. Nonlinear wave equation with nonequilibrium observables 



We now consider a nonlinear wave equation with a nonequilibrium set of observables. As will be seen, a surprising 
result is that the choice of observables forces us to nonlocally coarse grain. The nonlocality of the coarse graining has 
very interesting implications with regard to the resulting induced RG flow. The (driven) equations of motion that we 
are considering are 



dtn = 



9t0 = TT + ttlMl 



(20) 



where (f> and tt are real- valued flclds. The driving now includes generalized forces in addition to "true" forces. By 
expanding around equilibrium solutions of V^(/> = we find that for each real-space position x. 



B = 



ai 
a2 



and 



/3i 
132 



(21) 
(22) 



This driving allows for more states in (0, 7r)-phase space to be accessible compared to the driving in equation (jlSp . 
This, in combination with the form of y, ensures that both and 7r-dependent observables are being considered. By 
using exponential discounting, we find that the diagonal operator of HSV's is given by (see Appendix E]) 



x{Pl\k\ 



4a 



a5|k| 



al\k\) 



1/2 



(23) 
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Wk does not have the HSV's ordered from largest to smallest, so it is not truly in balanced coordinates. It is immedi- 
ately apparent that the HSV's are large for both large and small magnitude wave- vector. A heuristic explanation for 
this strange result is that for large wave-vector, tt is a pathologically "fast" variable. However, by driving tt with ui 
over all wave-vector, this permits the fast resonances to be excited at large wave-vector. The pathological nature of 
TT as an observable is analogous to the pathological nature of considering ^ an observable where ^ satisfies a Langevin 
equation [30]. In this sense, tt is a nonequilibrium observable of sorts. Furthermore, because both the small and large 
wavelength physics contributes strongly to the response of the system, local coarse graining cannot be correct. The 
appropriate coarse graining is nonlocal. 



IV. RG ANALYSIS 



Now that we have specified how to coarse grain, we must specify the distribution of initial conditions that we wish 
to coarse-grain. Should we choose a Gaussian distribution of initial conditions, then it turns out that provided that 
the action, S (g, {(/^o}), does not have a non-local spatial dependence, then the integration regions decouple in Fourier 
space and the RG is trivial. We thus consider a distribution of initial conditions that is local and quartic. Specifically, 
we consider the action for the initial conditions to be 



5(g,{0o})- J dxQv</>oV0o 



It is important to remark here that we consider the gradient term to be local (or marginally local). Equivalently, in 
Fourier space, the action is given by 



4! 



. 4 / 4 \ 

T J n l,^''-'' ) '^o(^i)'^o(''2)0o(k3)<^o(k4), (24) 



where D is the spatial dimension of the system. Rather than being wholly unmotivated, with the addition of kinetic 
terms, the CTP-method [1, [2^ may be used to derive a dynamical action from the above S that is approximately the 
same as Sdyn in Equation pO|) . Also, in the remainder, to avoid unnecessary subscripts, we denote (po by (p with the 
understanding that we're performing the RG on the distribution of initial conditions. 



A. RG equations from local coarse graining 

Here we introduce a large wave- vector cutoff fl. We define = 0< -I- 0>, where (/i< is only nonzero for |k| < Q/b 
{b > l)and 0> is only nonzero for fl/b < |k| < ft. Integrating out wave-vector shells between il/b and entails 
integrating out the 0> fields. We then rescale by defining 

0<(k) -Z(^i(&k) (25) 

and p = &k. Although we do not start with a "mass" term in the action (i.e. rn^0^), such a term is generated by the 
RG flow. Upon following the Wilson RG variant of what is described in Appendix [B] or the analysis done by Shankar 
in 0, we obtain the RG equations 

= (4 - D)\ - + OCX') (26) 

dim^ = 2rn' + ^^A(l - m^) + 0(P). (27) 

A and are dimensionless and result from appropriately rescaling A and by f2. The RG flow and fixed points 
for these equations are well studied. Two fixed points are the Gaussian fixed-point, ffi^ = and A = 0, and the 
Wilson-Fisher fixed-point, sa —(4 — D)/6 and A w 2(4 — D){2tt)^ / {3Sd)- Following standard convention, we 
approximate the Wilson-Fisher fixed-point in powers of 4 — Z? (the e-expansion) . We provide the form of these RG 
equations because, as will soon be evident, they differ greatly from those that we derive in the next section. 
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B. RG equations from nonlocal coarse graining 



In the case where and tt are treated on equal footing as observables, which in general may not be the case, 
ai = a2 — a and Pi — = (3. In the remainder, we treat this particular case. Furthermore, without loss of generality, 
we set a = (3 = 1. In this case, equation (|23p indicates that the |k| = 1 states are the least important. Thus, for the 
purposes of the RG, the |k| = 1 hyper-surface serves as our analog of the wave- vector cutoff. Implementing the second 
step of the procedure for generalized RG involves integrating out k-shells away from the |k| = 1 surface. Rather than 
transform the system into the balanced ^-coordinates, out of convenience, we coarse grain the system in k-space. 

In order to coarse grain, we let (j) ~ 4>< + 4>m + 0> where (j)^ is only nonzero for |k| < A, is only nonzero for 
A < |k| < A~^, and (/)> is only nonzero for |k| > A^^, where A < 1. Using this decomposition, the path integral 
measure factors as Dcj) = T)(j)^D(j)m1)(j)>- The RG equations are then induced by integrating out 4>m and then rescaling 
the wave- vectors and fields. For this problem, the rescaling procedure requires that 



and p = A^^k for |k| < A and p = Ak for |k| > A~^. Naive power counting breaks down as a direct result of rescaling 
in the two disjoint wave-vector regimes. 

Although we start with a theory where g = (1, A, 0, . . .) we can expect that the RG transformations may generate 
new nonlinear terms and that the coupling constants may become coupling functions. In fact, g flows towards having 
an infinite number of nontrivial components. In particular, in a complete treatment, the coupling constant A becomes 
a coupling function, A(pi, P2, Pa, P4). However, if we focus purely on the constant contributions of A(pi, P2, Ps, P4) 
in the different wave- vector scaling regimes, we see that it may be decomposed into the couplings {A(j^4_j)}^^]^. Here 
A(j j) represents the coupling constant in the scaling-regime with i wave- vectors having |p„| < 1 and j wave- vectors 
having |p„| > 1. We denote the mass terms for |q| < 1 and |q| > 1 respectively by and m?,. If we let A = e"*^', 
then to first loop order the RG equations for {A(j 4„j-)}^^j are 



0<(k)-Z<(^i(A-ik), 
0>(k) = Z><^2(Ak), 



(28) 
(29) 




(30) 



(31) 




(32) 



(33) 




(34) 



The RG equations for the mass terms are 




(35) 



(36) 
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bl < 1 






bl > 1 






coarse 




^™,n > 0, 


r, 




^™,n > 0, 




graining 


n>2 


m > 4 


n > 4 


n>2 


m > 4 


n > 4 


local 


no 


no 


no 


no 


no 


no 


nonlocal 


no 


no 


no 


yes 


yes 


yes 



TABLE I: Relevance of perturbations when D=4 



The first thing to notice in equations is that the contributions from tree level, the linear terms, indicate 

that the couplings involving a mixing of wave- vectors (i.e. z, j ^ 0) are irrelevant. The terms that are relevant, appear 
similar to equations (|^ - (P7)) . In fact, if we expand about m\ = 0, make the identifications A(4_o) = ^ and "^^r = "^^j 
and set the rest of the coupling constants to zero, then Equations (PO)) and ([55)) are exactly the same as equations 
(PS)) - (I?7)) . Thus, if it were possible to ignore all other relevant perturbations, then the new set of RG equations would 
retain the standard fixed points (Gaussian and Wilson-Fisher). With local coarse graining, higher-order derivatives 
and nonlinearities are rendered irrelevant by coarse graining. This then justifies why we only need to renormalize 
the fh? and A couplings. However, upon non-locally coarse graining the system, it is no longer possible to ignore the 
wave-vector dependence that A acquires. Such perturbations are relevant for |p,;| > 1. Specifically, higher derivative 
perturbations, p"(/)^,n > 2 and p" > 0,r7i > 4, in addition to higher order nonlinearities, (/)",n > 4, become 
relevant when |pi| > 1. While the couplings at small wave- vector, < 1, almost obey the standard RG equations 
obtained by local coarse graining, if the system is perturbed to its large wave-vector regime, then the couplings 
will flow away from the Gaussian or Wilson-Fisher fixed points [2|. This reflects that the dynamically faster short- 
wavelength perturbations are able to excite the conjugate field tt, thereby driving the system away from its standard 
statistical equilibrium. Were the conjugate field not accessible to the "noise", ai = 0, or not an observable, /?2 = 0, 
this would not have occurred. We summarize the relevancy of perturbations in Table U in the case when D = A. A 
key observation to make from our analysis is that the nonlocal coarse graining produces RG equations that are very 
different from the canonical ones. Equations (P5)) - (P7)) . In fact, a complete analysis would require us to determine the 
RG flow of infinitely many coupling constants. 

While the existence of an infinity of relevant perturbations renders Equations ((30)) - ((36)) meaningless, for complete- 
ness, we present their fixed point structure anyhow. One fixed point of the equations exist where « — (4 — D)/6, 
A(4,o) ~ 2(4-D)(27r)^/(3S'r,), to> ~ (4-D)/6, A(o,4) ~ -2{A- D){2tt)" /(SSd), and the rest of the couplings zero. A 
problem with it is that it doesn't guarantee positivity of the action. This already suggests the need to include higher 
order, relevant terms. For instance, to fix the positivity issue, we would need to add at least a (j)^ term to the action. 
Similarly, for -D < 4, the irrelevancy of the rest of the terms ensures that we would lose positivity, thereby 

requiring the addition of higher order terms. However, for D > 4, there is an alternative fixed-point that guarantees 
positivity of the action. At this fixed point, the coupHngs are m?, « (4 — -D)/6, A(o,4) « —2(4 — _D)(27r)^/(35_D), 
and the rest zero. While this fixed point seems strikingly similar to the Wilson-Fisher fixed point, it is important to 
recognize that the coupling constant multiplies products of fields in a different wave-vector regime. 

V. CONCLUSION 

In this paper we have presented a new RG procedure and have applied it to a 0^ toy model. We have shown 
that when both equilibrium and nonequilibrium observables are chosen, this RG procedure predicts that naive power 
counting breaks down and that terms that are ordinarily irrelevant become relevant at large wave-vector. We have 
shown that the RG equations that we have derived using non-local coarse graining differ significantly from the RG 
equations derived local coarse graining. While our RG equations superficially retain the Gaussian and Wilson-Fisher 
fixed-points as solutions, the infinity of relevant perturbations ensure that such fixed points either do not exist or that 
they pick up an infinity of unstable directions in the RG flow. Additionally, the equations superficially admit many 
other fixed points. The generalized Wilsonian RG developed here is applicable to nonequilibrium and heterogeneous 
systems, finite or infinite dimensional systems, and systems with various perturbations and uncertainties. Although 
the RG is still formally an uncontrolled approximation, the coarse graining is chosen such that the effective, coarsened 
system is close to the original one. Despite the versatility of this method, it is often difficult to analytically determine 
the balancing transformations. However, since there are very efficient numerical algorithms for finding balanced 
coordinates, this generalized RG remains a numerically useful and practical algorithm. 
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APPENDIX A: BALANCING X AND Y 



In this appendix we intend to calculate the balanced form of the operators X and Y (from Equations ([6]) and ([7])), 
called gramians, for linear wave equations (or their discretizations). This entails calculating the damped (exponentially 
discounted) gramians: 



As in the body of the paper, 



y(a) _ 



(Al) 



B = 



ail 
a-il 



Pil 
/32l 



Let us first introduce the following notations and conventions. Recall that any matrix, S, may be expressed in 
terms of the canonical matrix units, e^. In other words. 



where each Sij is just a complex number. For instance, in the case of 2 x 2 matrices. 



ei2 



1 




Additionally, for this section, Q 
instance, suppose 



then 



First note that if we define 



then easily it follows that 



1 
-1 



Lastly, we frequently use the algebraic tensor (dyadic) product, ®. For 



A- 



A®B 



All Ai2 
A21 A22 



AiiB A12B 
A21B A22B 



R = en ® _^ ^ ^1/2 



(A2) 



A = ei2 (g) I - 621 «) fl^ 



^ R^AR = Q ® rj 



I 



D. 

~n 



(A3) 
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From this, one then finds that 



R 



ajn 

ajfl-^ 



e-^®"*dtR 



However, using that e**^^^* = I ® cos Jit + Q (g) sinfl we finaUy arrive at 
X^"^ = R 



afJlcos^ nt + a^n^^ sin^ nt ^n^^iajn-^ - a?f2)^ sin^ nt 



dt R 



4a 









ajn + ajn- 



RH R X 

4a 



-a2n-l(a2l _ c.2j^2)(„2j ^ f22)-l ^(^2j _ a2f22)(„2j ^ ^^2)-! 

a(a2l ~ aln^)(an + n^)'^ a^n~\all ~ aln^)ian + O^)-! 



R. 



Similarly for the other gramian we obtain 



4a 



/3^n + pln-^ 

/3|f2 + /J^n-i 



R + ^R 

4a 



a{(3fl - Pin^){an + -a^n-\f3fl - 0in''){a^\ + O^)-! 



R 1. 



(A4) 



(A5) 



(A6) 



From equations (jASP and IA6I it follows after using 11^ to diagonalize $7 and taking the small "a" limit that the 
balanced gramian, without ordered eigenvalues, is given by 



4a 



1 







1/2 



/3|Ao + /3?A^^ 





/?|Ao+/32A, 



1/2 



APPENDIX B: RG EQUATIONS TO 1-LOOP 

In this Appendix, we briefly derive the RG equations for the two-point (Eqs. ([35]) and ([36])), and the four-point 
coupHng functions (Eq. ([30]) ). In Section [BT] we sketch out the standard perturbative procedure [U, [3, 3 used 
to derive the RG equations. Then we proceed to use the perturbative procedure to derive the RG equations for the 
two-point and four-point couplings in Sections IB 21 and IB 31 respectively. 



1. Sketch of perturbation theory 

Recall that we will be integrating out the field variables and keeping (/)< and 0>. Ideally we would like to 
exactly evaluate the partial trace 



2)0™ exp I -S'(0< , 0> , ) I = exp I -S'(0< , 0> , ) I 



(Bl) 



In general, computing such a functional integral is very difhcult. Consequently, we will perturbatively evaluate 
Equation (|B1[) about the quadratic part of the action. Now let us denote the quadratic part of the action by Si and 
the remainder by Sr (i.e. 5 = 5*2 -t- Sr)- By the simple form of the quadratic part of the action in Equation (|24|) . we 
have that 



s. 



4>>i4>m) = S2i4><,4>>) + S2{,(i>m)- 



(B2) 
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Now given that Vo{4'm} = m ^'^P{~'^2(<^m)} is the probability distribution generated by S2{<j>m), let (. . .)o denote 
averages taken with respect to this distribution. Using this notation, Equation (|Bip becomes 

T)4>„i cxp I -5((^< , , 0„O I = exp I -5*2 ((^< , 0> ) I (exp {-S^Dq . (B3) 

Note that we have absorbed any contribution from Zg^m into the path integral measure. 

Equation (jB3p may be (approximately) calculated by expansion techniques. Two techniques that are frequently 
used are the cummulant expansion T, '3| and Feynman diagram expansions [4| . Although Feynman diagram expansion 
methods permit infinite resummations more easily than the cummulant expansion, we make use of the cummulant 
expansion. However, we will still illustrate the non-vanishing Feynman diagrams that contribute to the cummulant 
expansion. Upon applying the cummulant expansion to Equation (|B3P we find that 



(expl-^r});) = exp 



(B4) 



The first term in the exponent on the right hand side of Equation (|B4p is the first cummulant while the second term 
(in square brackets) is the second cummulant. If we take S' to be the exponent on the right hand side of Equation 
(|B4I) then the Feynman diagram representation of S' is given by 



S' = 



6-6 




(B5) 



where the solid propagators represent the 0< and 0> propagators while the dashed propagator represent the (pm 
propagator. We will use the first and second cummulants (and diagrammatics) to derive the RG equations in the 
following sections. 



2. Two-point couplings 



If we consider having a mass term (i.e. — > — m^), then after integrating out 
according Equation the bare two-point propagator (Green's function) becomes 



and rescaling 0> and 



(B6) 



To simplify notation we denote the free propagator by Go{p)- Similarly, for \p\ < 1 we denote the propagator by 
Gq{p) and we do similarly for \p\ > 1. At one loop we obtain 



IT (a ' [A(4,o) /a<|p|<i GoiP) + \2a) k<\p\<x-^ (p)] 4|<i 1^(9)1' 
+A2 [A(2,2) /a<|,|<i G<{p) + A(04) /i<|p|<,-i G>{p)] |^(<z)p) 



/m<i 1^(9)1' + 



^(2,2) 



(0.4) 
1 



4|>il<^(9)l 



(B7) 
(B8) 



where Goip) is the free propagator and So is the area of the D-dimensional unit sphere. 
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and ([M]) ) 



By combining Equations (jB4p . (|B5p . (|B6p . and (|B7|) and using A = e we arrive at the RG equations (Eqs. ([35 

am^ - 2m2 + ^(4-°) + A 



V2,2) 



(9;m?, — —2m?, + 



2(27r) 



\(0,4) 



0(A2) 
0(A2) 



(B9) 
(BIO) 



3. Four-point couplings 



The decoupUng of the wave- vector regimes, while explicit in the calculation of the 2-point couplings, is more subtle 
in the calculation of the 4-point couplings. However, it is still manifestly evident at tree level. Consequently, as in 
the previous section, we will first calculate the rescalings for the bare 4-point vertex. 



X 



4 '^(4.0) 



4!(27r)'iD i|p,|<l Hi X2¥p-(27r)^ 



(Bll) 



X 



X 



\2(D-1) ^{3.1) f T-r3 f dPp±(OTr\D 

'k\{2Tr)io J|p.|<illi=l (2^)15 J|p4|>l (2-k)0\^^) 



^ E, = lPi V'iPl)viP2)viP3)v{P4:), 



(B12) 



= A^ 



4\{27r)*o J\pi\<l 



T-r2 rf°p. r 
iU=l (27r)0 J|pfc|>l 



d Pk 
3 (2Tr)0 



(27r) 



^ = ) </'(Pl)'^(P2)</3(p3)'^(P4), 



(B13) 



X 



X 



\2 -^(1.3) r T-r4 rf°Pi r d°Pi /o„W 

4!(27r)-4" J|p,|>lll4=2 (2ir)" J|pi|<l (27I-)" V^'^-' 

(Ej = lPi) <P(Pl)'^(P2)¥'(P3)¥'(P4), 



_ A4-D ^(0.4) r TT d^Pi C2t'I^ 

^ 4!(27r)'io J|p,|>l lli (27r)0 V^'^'' 



(B14) 



(Ej = iPj) 'y5(Pl)'^b2)</5(P3)¥'(P4). 

At one loop, a sample calculation with only 0< on the external legs of the 4-point diagram produces 



(B15) 



[^f (2^)-4^ /a<|q.|<a-i ^(91' 92, Qi,Q2)i27r)^6^ (gi + 92 + Qi + Q2) 

Afe, 94, Qa, Q4)(2^)^<5^ (93 + 94 + Q3 + Q4) (27r)^5^ (Qi + Q3) Go(Qi) 



(4! (2 J; 



(27r)^<5^ (g2 + Qi) GoiQ2MqiMq2Mq3Mqi) 
2 



(B16) 



)0/2 



/l 



<?.|<1 



/a<|Qi|,|Q2|<1 -^(4,0)^0 (<9i)Go(Q2) 
+ 2/a<|Qi|<i/i<|Q2|<A-i ^(34)^0 (<3i)G5"(Q2) 

'^■^i<IQiMQ2l<i •^(2,2)'^o'(*3i)Go (*32) 
{qi +q2+Ql+ Q2) {2tt)^S^ (93 + 94 - Ql - Q2) (y5(qi)^(g2)(^(g3)<y5('74). 



(B17) 
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From equation (|B17p . there is an acquired wave- vector dependence. For the sake of comparison with equations, we 
only wish to derive the (naive) flow equations for the constant part of X{qi, q2, Qs, 94)- Thus, we obtain 



Soil 



^(4,0) (l+"^<) 



+2X1^^^ {I + ml)-' {1 + ml)-' +Xl^^^ {l + ml)- 



(B18) 



Combining Equations (|B4l) . ([BS]) . (|Blip - (|B15P , and (jBlSp and using A = e yields the RG equations given in 
Equation ((50)) . 
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